BROKEN BRACELETS, MOLIEN SERIES, PARAFFIN WAX 
AND AN ELLIPTIC CURVE OF CONDUCTOR 48 



TEWODROS AMDEBERHAN, MAHiR BiLEN CAN, AND VICTOR H. MOLL 

Abstract. This paper introduces the concept of necklace binomial co- 
efficients motivated by the enumeration of a special type of sequences. 
Several properties of these coefficients are described, including a con- 
nection between their roots and an elliptic curve. Further links are given 
to a physical model from quantum mechanical supersymmetry as well as 
properties of alkane molecules in chemistry. 



1. Introduction 

A jeweler is asked to design a necklace consisting of a chain with n place- 
ments for k pieces of diamond. The client ask for one group of r diamonds 
to be placed next to each other and the remaining diamonds are to be iso- 
lated, that is, each one is mounted so that the two adjacent places are left 
empty. These special diamonds are called the medallion of the necklace. 
Figure [T] shows a necklace of length 20, with a medallion of length 5 and 
four extra diamonds. 




Figure 1 . A necklace with a medallion. 

Throughout this paper, a necklace is understood to be cyclically symmet- 
ric (unlabelled) formed by diamonds of two colors (i.e., binary). 
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Figure 2. A configuration. 



Figure 3. A forbidden configuration. 

A configuration or broken necklace is one resulting from one of the r + 
1 cuts to the left, right or in between the medallion (as described in the 
motivational introduction). Figure [2] shows a configuration and Figure [3] 
depicts a forbidden cut. 

Label n vertices as {1, 2, • • • ,n - 1, n}. The neighbors of the vertex i are 
i - 1 and i + 1 for 2 < i < n - 1 ; the single vertex 2 for i = 1 and the single 
vertex n— 1 for i = n. Configurations consist of a linear array of n vertices, k 
of which are marked or painted red. The marked vertices are either isolated, 
that is, its neighbors are not marked or connected, that is, the sequence of 
vertices [i, i + 1, i + 2, • • • , j} are all marked. In the latter case, it must be the 
case that i = 1 or j = n; that is, connected marked vertices contain 1 or n. 
Question 1. Determine the number Zk(ri) of configurations up to symmetry. 

The problem above, sans restriction, may be interpreted as a periodic 
chain made of two kinds of beads. The classical result on counting all 
necklaces with n beads is given by MacMahon formula 

(1.1) N(n) = -Y tp(d)2 n/d , 

where the summation runs through all divisors d of n, and <p(d) is the Euler 
totient function counting the numbers 1, 2, . . . , d relatively prime to d. 

Onofri et al [5] introduced a supersymmetric quantum mechanical model 
for a system whose degrees of freedom are bosonic and fermionic (creation 
and destruction) operator matrices. In the Hilbert space of this model, the 
vectors can be put in one-to-one correspondence with binary necklaces (or 
periodic linear sequences), where 1 and represent fermionic and bosonic 
matrices respectively. Pauli's exclusion principle provides a Fermi statistics 
which projects out a subset among all necklaces. As such the resulting space 
allows a genuine depiction of supersymmetry. 

The Hilbert space of the above model is generated by states due to sin- 
gle trace operators. Since trace is cyclically symmetric, all length n neck- 
laces related by cyclic shifts are identified as having the same state of n 
quanta. The Pauli principle is reflected in the fact that the fermionic op- 
erators are Grassmannian variables, hence they anti-commute; while their 
bosonic counterparts commute freely as scalars. The number of necklaces 
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is enumerated by MacMahon's formula < | 1 . 1 1 ) . Some of these are excluded 
as a result of anti-symmetry of planar states, dictated by supersymmetry. 
The following serves as an illustrative example with n = 4. There is a total 
of 6 necklaces listed here as linear sequences of period 4, namely 

0000, 0001, 0011, 0101, 0111, 1111. 

Supersymmetry breaks up this family into forbidden and allowed necklaces 
as follows. Given a sequence, start shifting a digit from right to left (prefix 
to suffix) and repeat until the original sequence is recovered. Every time 
a fermionic operator (that is, a 1) crosses another then a sign change must 
be registered because of anti-commuting; each bosonic operator (i.e., a 0) 
shifts around without any effect. At the end of the procedure, if the sequence 
becomes its own negative then call it forbidden, otherwise it is allowed. One 
of the main results of [5 J predicts 

(1.2) N allowed (n) = \Yj ^ T ' d ' and N forbMen(n) = \Yj V^ 2 '^ ' 

d\n d\n 

d odd d even 

Let us take a look at each of the above sequences now: 0000 stays the 
same; 0001 1000 0001; 0011 -1001 +1100 0011; 0101 
-1010 -0101; 0111 +1011 +1101 +1110 0111; 1111 -> 
-1111. Therefore there are four allowed {0000, 0001 , 001 1 , 01 1 1 } and two 
forbidden {0101, 1111} sets of necklaces. This agrees with 

Avowed = ^(1)2 4 ] = 4, and AWidden = ^(2)2 2 + <p(4)2 l ] = 2. 
The interested reader should find the complete story in J5]|. 

The central object of the work presented here is a sequence of numbers 
labeled necklace binomial coefficients Q . These coefficients have proper 



ties similar to the usual binomial coefficients: recurrences (Corollary 2.7), 



symmetries (Corollary 2.12) and an explicit formula in terms of the bino- 



mial coefficients (Theorem 2.8). Section [3] discusses a surprising result on 



the zeros of the necklace polynomial N t (y) (this is the generating function 
of )• It turns out that all its zeros are inside an elliptic curve. Moreover, 
this curve is the same for all values of t. Section |4]presents arithmetical and 
geometrical properties of these polynomials via generating function meth- 
ods. The necklace binomial coefficients also appear in a physical model of 
Onofri et al and in the study of symmetries of paraffin molecules in |@]|. 

2. The number of configurations 



In this section the counting problem from the Introduction is rephrased 
and solved. The current format as well as the original formulation will be 



4 



T. AMDEBERHAN, MAHIR BILEN CAN, AND VICTOR H. MOLL 



used interchangeably: 

determine the number Z k {n) of painting k points in red from a linear array 
ofn of them, with the condition that consecutive red points can only appear 
at the beginning or at end of the array. Moreover, arrays that are reflections 
of each other should be counted only once. 

In order to determine the number of configurations Z k (n) it is convenient 
to begin with a simpler count. 

Proposition 2.1. Let f k (n) be the number of arrangements of n vertices 
with k marked vertices, no consecutive marked ones where reflections are 
not identified. Then 

In - k+ 1^ 
(2.1) /*(*) = ( k 

Proof. Each such arrangement can be obtained by placing the k marked 
vertices and choosing k — 1 places to separate them. The count is obtained 
by eliminating the separating spaces. □ 

Reduced configurations. The next step is to count those configurations ob- 
tained by cutting the necklace exactly on one side of the medallion. These 
produce linear arrays where clustered vertices appear either at the beginning 
or at the end of the array. Invoking symmetry, only those with the medallion 
at the left will be considered. 

Definition 2.1. The function f3 k (n) denotes the number of linear arrays ob- 
tained by cutting a necklace with n vertices and k marked vertices with a 
medallion at the left of the array. 

Theorem [278] provides an expression for the function f3 k (n) and Theorem 



2.16 provides a formula for Z k (n). 



Definition 2.2. Let g k {n) be the number of arrangements of n vertices with 
k marked points, no two being consecutively marked and identifying sym- 
metric pairs. 

Example 2.2. A numerical example of g k (n) is given here. Take n = 4 and 
k = 2. From the pairs {12, 13, 14, 23, 24, 34} eliminate {12, 23, 34} for 
being consecutive somewhere. This leaves {13, 14, 24}. The pairs are now 
considered modulo 5, so that 24 is identified with 13 (the same as 31). The 
final allowed list is {13, 14} showing that gjfA) = 2. 

Theorem 2.3. The function /3 k (n) satisfies 

(2.2) p k (ri) = g k (n) + £ f k _ r (n - r - 1). 

r=2 
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Proof. Separate the different configurations into two groups: those with no 
consecutive marked points and those with at least two consecutive ones that 
are marked. The first type is counted by gk(n). Observe that if a certain 
arrangement has two or more adjacent marked vertices, then the remaining 
marked ones have no restrictions due to symmetry. In other words, reflec- 
tion only imposes limitations if the configuration has no adjacent marked 
vertices in it. 

The number of possible consecutive marked points is given by the size 
of the medallion. If this size is r, with 2 < r < k, then drop r + 1 places 
from the configuration (r for the medallion and one more at the right-end 
of it). This leaves a total of n - r - 1 spaces where to place k - r marked 
vertices. □ 

The next step is the enumeration of gk(n). This group is divided into 
three disjoint subclasses, those with (1) both ends are marked, (2) both ends 
are unmarked and (3) only the left end is marked. In the first class drop 
the vertices at positions 1, 2, n - 1 and n and observe that the remaining 
n - 4 vertices have k-2 marked ones and no further restrictions. Therefore 
there are gk-i(n -4) such arrangements. Similarly, the class (2) has gk(n - 2) 
elements. Finally, in class (3), drop the first two vertices and the last one that 
is not marked. The remaining n — 3 vertices have no symmetry restriction. 
The latter are counted by / fe _i(n-3) = ("^j 1 ) such arrangements. This gives 
the relation 

(2.3) g k (n) = g k {n - 2) + g k _ 2 (n - 4) " ~ k ~ ' 



k-\ 

Theorem 2.4. Let n = m + 2k - 1 and define gk(m) := gkiyn + 2^-1). Then 
gk satisfies 

(m + k — 2\ 
k- 1 



(2.4) 



gk(m) = gk-iim) + g k (m - 2) + 



Proof. Observe that any valid arrangement counted by gk(n) must satisfy 
n > 2k - 1 . The rest is elementary. □ 



The next result was obtained from experimental data generated by (2.4). 
Example 2.5. The function gk(m) is computed for < m < 3: 



(2.5) 
and 

(2.6) 



g k (0) = 1, g k (l) = 



k + 2 



, gk(2) = 



;=o L 



j + 2 



+ 



(k + 2) 2 
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Definition 2.6. The relation (2.4 1 attains a cleaner form by introducing the 

necklace binomial coefficients 



(2.7) 




for < k < t 
otherwise. 



The next result is a restatement of Theorem 12 .41 

Corollary 2.7. The necklace binomial coefficient satisfies the Pascal-type 
relation 



(2.8) 



k) n 



t-2 
k-2h 



+ 



t-2 
k-l 



+ 



t-2 
k 



The evaluation of the necklace binomial coefficients is now easy to guess 
and establish using ( |2.8[ ). 

Theorem 2.8. For < k < t, it holds that 



(2.9) 



Moreover, 



(2.10) 




for t even and k odd, 
otherwise. 



PkV) 



k+1 

k 



+ 



EC 



r=2 



Table [2] shows the values of the necklace coefficients: 
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8 


9 


10 
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1 














5 
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10 


9 
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12 


19 


19 


12 
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1 








8 


1 


4 


16 


28 


38 


28 


16 


4 


1 






9 


1 


5 


20 


44 


66 


66 


44 


20 


5 


1 




10 


1 


5 


25 


60 


110 


126 


110 


60 


25 


5 


1 



This array is known as the Losanitsch 's triangle and information about 
it can be found in Entry A034851 of Neil Sloane's Encyclopedia of Integer 
Sequences. 



A series of elementary consequences of (|2.9[) are presented next. 
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Corollary 2.9. The row- sum identity 

(2.11) V(') =2 f - 1 +2 L(? ~ 1)/2J 

holds. 

The next statements employ the Fibonacci numbers F n , defined by the 
relation F n = F„_i + F n _ 2 with initial conditions F = F { = I and the Lucas 
numbers L n defined by the same recurrence and with initial conditions L = 
2, L x = 1. 

Corollary 2.10. Let F„ and L„ as above. Denote f := [t/2j + 2 + 
Then 

t 



(2.12) Y/3 k (t)=UL t+2 + Fr)-l 

k=o 1 

Corollary 2.11. The generating functions 

(2.13) Yj^f = \{l+y)'+ l -{l + y 2 r 2 \l+y) 

(i i4) V (*) y = (i + ^) L(/:+1)/2J + (i-^ +1)/2J 



t mod 2 



2(l-x)^D/21(l-x 2 )L(' 



*+l)/2J 



and 



hold. 

Corollary 2.12. The necklace binomial coefficients are symmetric, that is, 

(2 - i6) CH-*. 

for < < 

Corollary 2.13. The function g is symmetric; that is, 
(2.17) g k {m) = g m (k). 



Proof. This is a restatement of (2.16). An alternative proof of the symme- 
try ( |2.17[ ) is obtained from the recurrence ( |2.4| ). Simply express it in two 
different forms 

(m + k — 2\ 

(2.18) g k (m) - g k - 2 (m) = g k (m-2) + 



<c(m) - g k (m - 2) = g k „ 2 (m) + 



k-l 
m-vk — 2 
k-l 
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The result now follows by induction and the symmetry of the binomial co- 
efficients. □ 

The next theorem provides a combinatorial proof of the symmetry rule 
4217T ). 

Theorem 2.14. The symmetry gk(m) = g m (k) holds. 

Proof. The assertion amounts to gk(m + 2k - 1) = g m (k + 2m - 1). Take a 
linear array of n nodes and its 2-coloring (red r or white w). By definition, 
gk(n) enumerates all possible ways of coloring k nodes in red with the rule: 
(1) no two reds are consecutive; (2) two such arrays are equivalent if they 
relate by reflection. According to (1), it must be that the first k - 1 reds are 
each followed by white. Thus, any selection of k reds can be interpreted as 
choosing the (k- 1) pairs rw and a free r. For each pair rw, trim-off the w as 
well as its sitting node. That means, when n = m + 2k - 1 then the number 
of nodes reduces to m + k and hence gk(m + 2k - 1) induces an equivalent 
counting of (m + fc)-nodes of which k are red (note: rule (1) is absent but 
rule (2) stays). Similarly, g m (k + 2m - 1) tantamount to the counting of 
(m + fc)-nodes of which m are white. But, it is obvious that coloring k nodes 
red on an (m + fc)-array is equivalent to the coloring of m nodes in white. 
This gives the required bijection. The proof is complete. □ 

Example 2.15. This example demonstrates the above proof; i.e. gk(m + 2k- 
1) = gm(k + 2m - 1). Take m = 2 and k = 3. Then, g3(7) and #2(6) count 
respectively the cardinality of sets 

A := {rwrwrww \rwrwwrw \rwrwwwr,rwwrwrw \rwwrwwr,wrwrwrw} 

and 

B : {rwrwww,rwwrww,wrwrww,rwwwrw,wrwwrw,rwwwwr}. 

The set B after color- swapping turns to 

Bi := {wrwrrr,wrrwrr,rwrwrr,wrrrwr,rwrrwr,wrrrrw}. 

The two sets A and B\ are now mapped (w-trimmed and r-trimmed, respec- 
tively) to 

A\ := {rrrww,rrwrw,rrwwr,rwrrw,rwrwr,wrrrw}, 

and 

Bn := {wwrrr,wrwrr,rrwrr,wrrwr,rwrwr,wrrrw}. 
The bijection between A\ and B u is clearly exposed; that is, reflect B u 
to get the set 

5m := {rrrww,rrwrw,rrwrr,rwrrw,rwrwr,wrrrw}. 
The full counting solution to the configuration problem is presented next. 
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Theorem 2.16. The total number Z^it) of possible linear configurations of 
k diamonds (with or without a medallion) on t nodes is given by 




Proof. Catalog the diamonds according to whether they are: (1) an equal 
number of clusters; (2) unequal number of clustered diamonds on the two 
end-nodes. However many are remaining to be mounted in the interior, case 
(1) is affected by the reflection but those in case (2) are not. It follows that 
the first case is enumerated by the function gk(t) (equivalently, by necklace 
binomials) while the function fk(i) is the right choice for the second cate- 
gory. The details are omitted. □ 



The necklace coefficients are given as Entry A005994 in Neil Sloane's 
Encyclopedia of Integer Sequences. The reader will find there information 
on the connection between \£} and the so-called paraffin numbers. The 
chemist S. M. Losanitsch studied in H the so-called alkane numbers (called 
here the necklace numbers) in his investigation of symmetries manifested 
by rows of paraffin (hydrocarbons). In the molecule of an alkane (also 
known as a paraffin), for n carbon atoms there are In + 2 hydrogen atoms 
(i.e. the form C n H 2n+ 2)- Each carbon atom C is linked to four other atoms 
(either C of H); each hydrogen atom is joined to one carbon atom. The 
figures in the Appendix show all possible alkane bonds for 1 < n < 5. 
There are 1, 1, 1, 2, 3 possible alignments, respectively. 

A geometric interpretation. Given a finite group G, it is a classical prob- 
lem to find the generators of the ring of polynomial invariants under the ac- 
tion of G. The Molien series M(z; G) is the generating function that counts 
the number of linearly independent homogeneous polynomials of a given 
total degree d that are invariants for G. It is given by 

1 1 °° 

(2.19) M(z; G) = — Y — = V b iZ \ 

\G\ £ det(/ - zg) £ 

Thus, the coefficients b\ record the number of linearly independent polyno- 
mials of total degree i. 



Now assume k = 2m - 1. Then (2.14) becomes 



vn/«' + 2m-l\ . 1 1 1 1 

< 2 - 20 > Z 2m- 1 / = 2(r#* 



i>0 v /9i 



2(l- z ) 2m 2{\-z 2 ) m 
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This is recognized as 

1 1 1 1 _ 1 y 1 

( — * 2 (1 - z) 2m + 2 (1 - z 2 ) m ~ \G\ j£ det(l 2m - zgY 

where G is the symmetric group S 2 and the summation runs through the 2m- 
dimensional group representation of the elements g in GL 2m (C). The argu- 
ment below shows that the series is indeed a Molien series for the ring of in- 
variants under the action of S 2 . More specifically, the ring of invariants un- 
der consideration is C[X; Y]~ 2 where X = (jci, . . . , x m ) and Y = (y u .. . ,y m ). 
The action is given by xi h-> ji for I = 1, . . . ,n. 

Let cr be the matrix cr = K ^ j and let n be the tensor product n = <t<S>1,„ 

resulting in a 2m x 2m matrix which has four blocks of size mx m with the 
off-diagonal blocks being the identity matrix and the diagonals blocks being 
zero. The matrix group generated by n in GL 2m is 5 2. Consequently, 

(2.22) det(l 2m - zn 2 ) = det(l 2m - zl 2m ) = (1 - z) 2m 

( 1 -z\ 

and det(l 2m - zn) = det(p <8> l 2m ), with P = \_ z j ■ Since det(A ® B) = 
det(A) m det(5) m , it must be that det(l 2m - zn) = (1 - z 2 ) m . 

These observations are summarized in the next statement. 

Theorem 2.17. Consider the action of Z 2 on C[x\ , • • • , x m ,y\, ■ ■ ■ , y m ] given 
by xi i-> yi. Then, the number of linearly independent invariant polynomials 
of total degree i is given by the necklace binomial coefficient (' 2 ^ 1 ) 3i - 



by 



3. THE NECKLACE POLYNOMIALS 

In this section we discuss properties of the necklace polynomials defined 



(3.D ™-zG)/ 

The explicit formula 
(3.2) N t (y) = i(l +yY + ~(1 +/) Lf/2J (l + y y mod2 



is given in (2.13 ) 
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Example 3.1. The first few values of N t (y) are given by 

MOO = l+y 

N 2 (y) = l+y + y 2 

N 3 (y) = N i (y)N 2 (y) 

N 4 (y) = 1 + 2y + Ay 2 + 2y 3 + y 4 

N 5 (y) = M(y)iV 4 (y) 

N 6 (y) = N 2 (y)(\+2y + 6y 2 + 2y 3 + y 4 ) 

N 7 (y) = N i (y)N 2 (y)(l+2y + 6y 2 + 2y 3 +y 4 ) 

N 8 (y) = 1 + Ay + \6y 2 + 28y 3 + 38/ + 28y 5 + 16y 6 + Ay 1 + y 8 . 

The sequence of necklace polynomials have some interesting divisibility 
properties. The results presented below began with the empirical observa- 
tion that, for t odd, N t (y) = Ni(t)N t -i(y). 

Corollary 3.2. Let j e N and t e N. Then Nj(y) divides N (2 t-i)j(y)- 

Proof. This is a direct consequence of the explicit formula given in Theo- 
rem 1X1 □ 

Problem 3.3. Prove that N^iy) is irreducible. 

Many polynomials appearing in combinatorics are unimodal; that is, there 
is an index n* such that the coefficients increase up to n* and decrease from 
that point on. A stronger property is that of logconcavity: the polynomial 
P(x) = 2Zk=o a k xk is logconcave if a 2 k - a k -ia k+ i >0forl < k < n - 1. The 
reader is referred to [0171 for surveys on these issues. 

The explicit expression (2.9) gives an elementary proof of the next state- 
ment. 

Theorem 3.4. The necklace binomial coefficients are unimodal. 
Proof. The inequality 

(33) {']<<• 



k L \k + 1 



for < k < lt/2\ and the symmetry of the necklace binomial coefficients, 
established in Theorem 2.12[ give the result. □ 

Theorem 3.5. The polynomial N t (y) is logconcave. 



Proof. Use (2.9 ) and separate cases according to the parity of t and k. □ 

Problem 3.6. Let 2{a n } := { } be an operator defined on non- 

negative sequences. Therefore, a polynomial P{x) is logconcave if £ maps 
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its coefficients into a nonnegative sequence. The polynomial P is called k- 
logconcave if 2^\P) is nonnegative for < j < k. A sequence is called 
infinitely logconcave if it is fc-logconcave for every e N. 

A recent result of P. Branden Q~| proves that if a polynomial P has only 
real and negative zeros, then the sequence of its coefficients is infinitely 
logconcave. The sequence of binomial coefficients satisfies this property. 

The question proposed here is to prove that N t (y) is infinitely logconcave. 

There is a well-established connection between unimodality questions 
and the location of the zeros of a polynomial. For example, a polynomial 
with all its zeros real and negative is logconcave [8J. This motivated the 
computation of the zeros of N t (y). Figure [4] shows the zeros of MooOO- 



-1.0 



-0.8 



-0.6 



-0.4 



Figure 4. The zeros of the necklace polynomial Niooiy)- 



Theorem 3.7. Let y = a + ibbea root of the necklace polynomial N t (y) = 0. 
For a ^ -1, define the new coordinates u = 1/(1 + a) and v = b/(l + a). 
Then (w, v) is on the elliptic curve v 2 = u 3 - 2u 2 + 2u - 1. 



(3.4) 



a+ y y = 



Proof. Any zero of N t (y) satisfies 

\(i+y 2 ) t/2 

[(1 +v 2 ) ( '~ 1)/2 (l +y) 

Taking the complex modulus produces |1 + y\ 4 
y = a + ib this equation becomes 

a(a 2 + a + 1) 



(3.5) 



b l = 



1 + a 



if f is even 
if t is odd. 

= |1 + y 2 \ 2 . In terms of 
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The transformation 1 + a = l/u and b = v/u leads to equation 

(3.6) v 2 = u 3 - 2u 2 + 2u - 1 =(u- l)(u 2 - u + 1), 

as claimed. □ 



Note 3.8. The collection of points on an elliptic curve (£, such as (3.6), has 
been the subject of research since the 18th century. The general equation of 
such a curve is written as 

(3.7) y 2 + a\y = x 3 + a 2 x + a^x + 

and if x,y e P(C 2 ), the complex projective space, then S is a torus. The 
addition of this torus is expressed on the cubic in a geometric form: to add 
Pi and P 2 , form the line joining them and define P 3 := Pi®P 2 as the reflec- 
tion of the third point of intersection of this line with the cubic curve. This 
addition rule is expressed in coordinate form: the general formula given in 
(6). Let P\ = (xi,yi) and P 2 = (x 2 ,y 2 ). Define 

{ y2-yi \f v -h y [ yi*a-yz-*i if v ■+ v 

— — It X 2 t Xi It X 2 ± X\ 

£?-4x 1+ 2 _ and v = \_*f +2 %_ 2 

— 2^ — itx 2 -x u [ — ^ — rtx 2 -x 1 . 

The P 3 = (x 3 ,y 3 ) is given by 

x 3 = A 2 + 2 - x\ - x 2 and y 3 = -Ax 3 - v. 

Aside from the point Pq = (1,0), the table below shows a collection of 
points on the curve (£ obtained using Mathematica. The notation 

y = >/3 + 2V3, 6 = yf/^2, t = >/24+ 14 V3, a = 2 ^2(11 +5 V5), 

&>i = 2 + V3, o> 2 = 2(3 + V5), a> 3 = 3 + 2 V3 

is employed. 

The notation necklace point refers to a point (u, v) on the elliptic curve 
($ that is produced by the zero y = a + ib of a necklace polynomial via the 
transformation 1 + a = l/u and b = v/u. The addition of two necklace 
points sometimes yields another one. For instance, P x ®P{ = P and 2P 3 : = 
P 3 © P 3 = P 2 . On the other hand, the set of necklace points is not closed 
under addition: 

p l& p 7 = I/7 + 3V5 + V 66 + 30 V5 j-^ ^21 + 9 V5 + -^30(29 + 13 V5) J 

The minimal polynomial for this number is y 8 - 28y 7 + 1948y 6 - 5236y 5 + 
4858y 4 - 3988y 3 + 7156y 2 - 6040y + 2245. This polynomial does not divide 
a N t (y) for 1 < t < 1000. It is conjectured that it never does. 
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Name 


u 


V 


Root of N t (y) = 


p. 


2 


-V3 


2 


r 2 


2 


+ V3 

~ V -J 


2 


Pi 


0)1 -y 


0)3 - T 


6 


Pa 


o)i - y 


-0) 3 + T 


6 


Ps 


0)1 +y 


-0)3 - T 


6 


Pe 


0)1 + y 


0)3 + T 


6 


Pi 


(1 + 8)0)2 


0)2 + CT 


4 


P S 


(1 + 6)o) 2 


-(cc>2 + cr) 


4 


P9 


(1 - 6)0) 2 


0)2~ CT 


4 


P10 


(1 - 5)0) 2 


-(0) 2 ~ Cr) 


4 



Table 1 . Some points on the elliptic curve (£. 



Note 3.9. Equation ( |3.5 ) shows that any root of N t (y) must satisfy -1 < 



Rey < 0. Observe that y = is never a root. 

Note 3.10. The change of variables « h « + 1 transforms the curve <£ into 
the form v 2 = u 3 + u 2 + u. This curve appears as 48a4 in Cremona's table 
of elliptic curves, available at 

|http : //www .ma .utexas . edu/users/tornaria/cnt/cremona .html? 



conductor=48 

The discriminant of the cubic is negative. Therefore the curve has a single 
real component. This is seen in Figure [4] 

Problem 3.11. The zeros of the polynomial N t (y) are algebraic numbers 
lying on the elliptic curve ([33]). The points on that curve with algebraic 



coordinates form a subgroup J?l under the addition described above. The 
question is to characterize in J{ the set coming from necklace points. 

4. Necklaces and their progeny 

This section explores the enumeration of certain special necklaces and 
their generating functions. The latter is applied to the computation of some 
Molien series. A circuit graph is a graph consisting of n vertices placed on 
a circle with some of them colored by red. 

Proposition 4.1. The total number of n-bead (circular) necklaces on which 
a red-red string is forbidden is given by 

(4.D "M = ;E*(5)** 

d\n 

Proof. A standard application of Burnside's lemma. □ 
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Example 4.2. For n = p prime, formula ( 4.1| ) gives 

(p-l) + L p 



(4.2) 



W(p) = 



p 2 W(p 2 ) = L p2 +(p- 1)L„ + p(p - 1). 



It follows that L p = 1 mod p. Similarly, for n = p 2 , (4.1 ) gives 
(4.3) 
It follows that 

(4.4) L p 2 = L p + 1 mod p 2 . 

These are well-known results (3). 

A more distinguishing count is provided by defining Wk(n) to be the num- 
ber of n-bead necklaces on which a red-red string is forbidden, consisting 
of exactly k red beads. In order to accomodate the possibility that k = 0, we 
define W (n) := 1 (this is justifiable since W (n) = \ X^^d) = 1) • 

Theorem 4.3. The function Wkin) is given by 

(4.5) ■ 1 



W k (n) = tYM** 

n - k 



d\n,k 



n k 

d d 
k 
d 



Proof. It follows directly from Burnside's lemma. 
Corollary 4.4. The identity 

L«/2J . ,n _ k 

d 

holds. 



□ 



k=0 d\n,k x d ' d\n 



Proof. The assertion follows from the combinatorial identity 
(4.7) Yj W ^ = 



k>0 



Theorem 4.5. For neN define 



□ 



(4.8) 



V d (x) = 



l - VTT4^ ( l + VT+l^ 



+ 



Then the row-sum generating function of Wk(n) is given by 
(4.9) F„(x) := J] W k (n)^ = - £ ^ V, . 



/t=0 
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Proof. The proof is based on the identity 

|m/2J 



1 L ^4 J 1 lm-k\ 

(4.10) -V m (x) = V 

m f— ' m — k\ k j 

k=0 v ' 

which is easy to verify. This is applied to 

L»/2J 1 / 

2><"^ - 2>>Z^ 



d \ v dk 
k 1 



k=0 d\n k>0 

Cf(d) ^ 1 (\-k 



Z <PW y _ 



, d f-i'j-k\ k 

d\n k>0 d V 

The result follows from here. □ 



Example 4.6. For p prime, the polynomial F p (x), defined in (j479j), is given 
by 

LW2J , I , \ 

(p - 1)2^(1 - Vl + 4jc) p + (1 + Vl + 4*y 
p-2P 



Example 4.7. Put n = 3k + 1 in (4.3 1 to obtain W k (3k + 1) = 2^l( 2 \ +1 ), the 
Catalan numbers. 

Example 4.8. For neN, and with L„ denoting the Lucas number, 

1 n-k\ 1 



(4.11) 



77 • 



This is obtained from setting x = 1 in (4.9). 



Theorem 4.9. The ordinary generating function for the diagonals of Wk(n) 
is given by 

„ 1 <P(d)x 2k 



(4.12) Yj Wk{n)xn = kYj 



n>k d\k v ' 

In its lowest terms, the denominator of this rational function takes the form 

(4.13) \\( 1 - x d f {kld) = Y] ®d(xf d , 

d\k d\k 

where O c/ (x) is the d-th cyclotomic polynomial given in terms of the Mobius 
//-function as ® d (x) = Il c |rf(l - x? lc Y (x) - 
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Proof. The result follows from the Taylor series expansion 

x2m v 1 (J- m \ i 

(4 - 14) "71 w = X " F- 

m( 1 - jc) m ^ j — m\ m ) 



A geometric interpretation. The above generating function W k (n)x n 
is the Molien series W(x;Z k ) for the ring of invariants C[X] Z * where X = 
(xi, . . . , X*;)- In this case, the group Z k is identified with its ^-dimensional 
group representation in GL k (C). More concretely, Z k = (e k ) where e k is the 
kxk permutation matrix such that e[i, j] = 1 if j = i + l;e[k, 1] = 1 and 
7] = 0, otherwise. Let RP(d) be the set of positive integers less than d 
and relatively prime to d. Partition the integer interval [k] into the disjoint 
union 

(4.15) [k] = {l,2, ...,k} = U^RP(d). 

d\k 

This relation is reminiscent of the well-known identity k = Y,d\k <p(d)- Then, 



1 k 1 

W(x;Z k ) = — V - 

|Z*l^det(l*-xe£ 



= -y 

* V det(l, - 
ly 



k det((l rf - xe rf ) ® \ kjd ) 



k^<mi d -xe d y<'* 



<p(d) 

These findings are stated in the next result. 

Proposition 4.10. The number of linearly independent homogeneous poly- 
nomials, of total degree n, for the ring of invariants C[X] Z * equals 

1 ^ /n + k 

1 . ~l d 



n + k Ti I § 

d\n,k x d 
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5. A SAMPLE OF THE COMPUTATION OF ZEROS 

Motivated by the interesting properties of the zeros of necklace polyno- 
mials, this section presents some computational graphics showing the ze- 
ros of the polynomials F n {x). Figure [5] shows the location of the roots of 
F ioooW- 



3 - 
2 - 




• • • ' -2 - 

-3 - 

Figure 5. The zeros of FioooC*)- 
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The next four figures show a selection of regions from the set of the roots 
of all the polynomials F n (x) for 3 < n < 1000. The caption indicates the 
range depicted. 




Figure6. [-10,2] x [-2,2]. Figure 7. [-2,0.2] x [-0.6,0.6]. 




Figure 8. [-1.4,-0.6] x [-0.5,0.5]. Figure 9. [-50,5] x [-6,6]. 

The interesting structure depicted in figures [6] to [9] will be explored in 
future work. 
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Appendix A. Rows of paraffin 

The figures show all possible alkane bonds (paraffin) C n H ln +2 for n = 
1, 2, 3, 4, 5. 
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